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ABSTRACT 

We present here a nonlinear Alfven soliton solution of the MHD 
equations. This solution may represent the final state of modula- 
tionally unstable Alfven waves. We also briefly describe a model of 
the expected turbulent spectrum due to a collection of such solitons. 


^Currently a Nuffield Foundation Science Research Fellow 
1 . Introduction 


The Alfvenic nature of the large amplitude turbulence in the 
solar wind has been known for some time (Belcher and Davis 1971). The 
relationship between this turbulence and traditional linear MHD waves 
or special large amplitude solutions (Barnes and Hollweg 1974; 
Abraham-Shrauner and Feldman 1977) remains unclear. Equally unclear 
are the consequences this turbulence may have for solar wind particle 
populations, although ion drift speeds proportional to the Alfven 
speed and high perpendicular temperatures suggest some wave-particle 
connections (Schwartz et al. 1981). 

Efforts to describe solau: wind turbulence usually begin with a 
perturbation expansion, leading either to the "decay cascade" (Chin 
and Wentzel 1972; Cohen and Dewar 1974; Schwartz 1977) if the system 
is described purely in terms of its linear modes, or to the 
modulational instability (Derby 1978; Goldstein 1978) of a single 
large amplitude Alfven wave. We discuss here a third alternative 
which leads to a nonlinear soliton solution. This solution may 
represent the state to which a modulationally unstable Alfven wave 
relaxes. Using an ensemble average of collections of these solitons, 
we find a turbulence power spectrum €(k)«k“2 at high frequencies 
which flattens toward the low frequency end. Details of the calcula- 
tions can be found in Ovenden et al. (1983). Further work on 
alternative turbulence models and particle— soliton interactions is 
currently in progress. 

2. The Soliton Solution 

We begin with the one fluid description of a plasma: continuity 
and momentum equations and the generalized Ohm's law coupled with 
Maxwell's equations (eg. Boyd and Sanderson 1969 3.40, 3.54 and 3.61— 
3.63). Assuming quasi-neutrality, an isothermal equation of state, 
and neglecting the displacen«nt current, resistivity and the 9j/dt 
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term in the generalized Ohm's law, it is eetsy to eliminate the 
electric field and transverse fluid velocity. For simplicity we 
consider only peurallel-propagation so that 3/0x - o - d/dy where 
the z direction corresponds to that of the background magnetic field 
Bq. Ultimately, the resulting equations can be written 
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®± “ Bjj ± iBy, V is the z— component of the fluid velocity, n is 
the density, vj^-Bq/ V( 47rp ) is the Alfven speed, fli the proton 

J^ 6 < 3 uency , Cg the isothermal sound speed, and d/dt™ 3/3t + v3/3z 

Both B^ and n have been normalized to their background values (Bq and 

n© respectively). The harmonic solutions to (1.1) and (1.2) with lB±l 
independent of z are just the Alfven— Ion cyclotron modes investigated 
by Abraham-Shrauner and Feldman (1977). 


We are interested in modes with small but nonzero fluctuations On 
and V in the density and velocity respectively. To this end, we 
linearize (1.1) and (1.2) with respect to these quantities. Further- 
more, we insert an Alfven wave with slowly varying amplitude, viz.. 
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where Qq* - cja ( 1 T cajv/zni ) and assuming u)i^<<n± with (Ja “ ^'a Ya* At 

this point, the modulational insted^ility calculation proceeds by 
perturbing the b - constant solution with fluctuations Ob, On, v all 
varying harmonically in space and time. Such a calculation yields 
instabilities with e-folding lengths “ 10”^ AU. 

Here, following Zedcharov (1972) we seeurch instead for soliton 
solutions of (l.l)-(i.2) with a wave amplitude in (1.3) given by 

b(z,t) - b^ sech [ x(z-Vt) ] ^ , (1.4) 

keeping terms to third order in I b I emd assuming b varies slowly 
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This process yields self- 
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for values of the ratio of sound to Alfven speed squared ( 0 ) << 1 and 
for the Wo+ ( 1®^* hand polarized ) wave . The density ( in non- 

dimensional units) and velocity vcuriations are given by 

lb(z,t) 1^ 
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Physically, the increased density (1.6) leads to a lower local 
Alfven speed causing the associated build up in wave energy in this 
region. The energy transport is achieved by forcing this material to 
move at the wave group velocity (l.5a), thereby resulting in the 
unchanged soliton envelope. It is not yet clear whether the solution 
(1.4)-(1.9) represents a true soliton, in the sense that two such 
solutions would pass through one emother without interacting. 


2. Turbulence Spectrum 


To build a picture of a turbulent spectrum out of solitons, we 

follow the work of Kingsep et al. (1973) and Yu and Spatschek (1976) 

who consider an ensemble of systems with a common energy density W 
(normalized to BQ^/en) and length L. The various realizations 
correspond to different numbers, N, of non-overlapping, assumed 

identical, solitons into which this energy is divided. The value of N 
lies in the interval 

1 < N < KL 

Ensena>le averaging over these realizations, which we take with equal 

probablities, yields an energy spectrum 

< lb. 1^ > « — - — f ak tanh (ok) - log cosh (ok) ] (2.1) 
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As can be seen in 
the plot of (2.1) shown 
in Fig. 1 (or directly 
from the explicit form of 
(2.1)), this spectrum is 
flat at low frequencies 
EUid falls off steeply as 
k“2 at high frequencies. 

Although this is 
qualitatively similar to 
solar wind observations 
(eg. Bavassano et al. 

1982 ; Denskat and 
Neubauer 1982), it is not 
yet clear whether 
quantitative agreement 
can be improved within 
the soliton description 
by, eg . , an alternative 
ensemble average, 
collections of unequal 
solitons , etc . 

3 . Summary 

We have found a nonlinear Alfven soliton solution which we 
envisage as the final state of a modulationally unstable Alfven wave. 
Dividing turbulent energy in a given region into such solitons and 
averaging over the different possible arrangements yields a turbulent 
spectrum which is flat at low k and falls as k~2 at large k. Remaining 
questions concern other possible descriptions of solar wind turbu- 
lence within this context and the consequences of particle-soliton 

interactions for solar wind ion distributions. 
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